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Abstract (Basic) : EP 624847 A 

The generator produces an approximation function based on pairs of 
values (Xi,Yi) associating a dependent magnitude (Yi) with an 
independent magnitude (Xi) , and generates new data pairs (Xa,Ya) 
following the approximation. 

The generator has a first part (10) that iteratively determines a 
current linear regression function to select those of the current 
regression fiinctions that deliver the approximation to the current set 
of input values with least error. The parameters (p,q) of the best 
regression function are then retained, and used by the second part (17) 
to generate the new approximate pair of data values. 

ADVANTAGE - Requires little computational capacity and gives rapid 
results under controlled maximum error, for use in neural networks or 
in approximation to measured values in industrial processes . 
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Abstract (Equivalent) : US 5519647 A 

An apparatus (5) for generating an approximation function based on 
first pairs {(XI, Yl) to (X6, Y6)) of values associating a dependent 
variable (Yl to Y6) with an independent variable (XI to X6), and for 
determining second pairs (XA, Y'A) of values of said variables in 
accordance with said approximation function, characterized in that the 
apparatus comprises: 

first means (10) : 

for iteratively determining at least one current linear regression 
function by making first errors of alternating sign equal in absolute 
value, which first errors have been measured between, respectively, 
first values (Y3, Y4, Y5) of the dependent variable for three pairs 
(X3, Y3) (X4, Y4) (X5, Y5) of a series of pairs, and second values 
(Y'3, Y'4, Y'5) of the dependent variable determined, in accordance 
with said current linear regression function, for the same values (X3, 
X4, X5) of the independent variable, 

for selecting that one of the current linear regression functions 
which produces an approximation of all the pairs of said series with a 
minimum error, and 

for coding the selected linear regression function with the aid of 
specific codes (p, q) , and 

second means (17) for determining said second pairs (XA, Y'A) with 
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© Dispositif (5) pour generer une fonction d'approximation fondee sur des prenniers couples ({X-.. Yi) a (Xs. 
Ys)) de valeurs associant une grandeur dependants (Y, a Ye) a une grandeur independante (X, a Xs). et pour 
determiner des seconds couples (Xa. Y'a) de valeurs desdites grandeurs d'apres ladite fonction d'approximation. 
Le dispositif comprend : 

- des premiers moyens (10) pour determiner iterativemant au moins une fonction lineaire courante de 
regression, pour selectionner celle des fonctions lineaires courantes qui delivre I'approxtmation de tous les 
couples de ladite suite avec des erreurs minimales, et pour coder, a I'aide de codes specifiques (p. q),'la 
fonction lineaire de regression selectionnee, 

- et des seconds moyens (17) pour determiner lesdits seconds couples (Xa. Y'a) a I'aide desdits codes 
specifiques. 

Le dispositif peut etre utilise pour le calcul de valeurs approchees de fonctions mathematiques. par exemple 
dans un reseau de neurones, ou pour determiner une fonction de regression formant une approximation de 
valeurs experimentales de mesure. par exemple des mesures dispersees issues du suivi d'un processus 
industriel. 

L'invention concerne egalement une methode pour generer une fonction d'approximation. 
Application : calcul approche de fonctions mathematiques. calcul de fonctions de regression. 
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L'invention concerne un dispositif et une methode pour generer une fonction d'approximation. cetle-ci 
etant fondee sur des premiers couples de valeurs associant une grandeur dependante a une grandeur 
independante. et pour determiner des seconds couples de vaieurs desdites grandeurs a partir de lad.te 
fonction d'approximation. 

Un dispositif et une methode de ce type sont connus du brevet US-A- 3 789 203 qui decrit un 
generateur de fonctions operant une approximation par interpolation iterative. Ce dispositif est prevu pour 
etre utilise dans des applications de traitement de donnees necessitant un calcul de fonctions teiies que 
sin(x). tg(x) par exemple. Ce dispositif ne requiert qu'une capacite de stockage minimale de la part d'un 
dispositif utilisateur. A partir de deux points appartenant a une fonction a interpoler. la methode tout d'abord 
interpole la fonction par une droite reliant les deux points, puis fait une approximation des ecarts entre la 
droite et la fonction par des approximations polynomiales de degre croissant. Ensuite elle substitue aux 
points initiaux des points approximatifs pour reduire la longueur du segment reliant les points a trailer el 
enfin reitere les operations precedentes. 

Une telle methode necesslte des ressources importantes en moyens de calcul et ne peut etre mise en 
15 oeuvre qu'avec des calcuiateurs performants. 

Or. il existe des applications ou une teUe methode ne peut pas etre mise en oeuvre car elles ne 
disposent pas des ressources suffisantes. De plus, pour certaines applications on peut se satisfaire d'un 
calcul approche de la fonction pour des valeurs, en nombre limite, de la grandeur independante. 

II peut s'agir d'une fonction sigmoVde appliquee h des potentiels neuronaux delivres par au moins un 
neurone dans un reseau de neurones, li peut s'agir d'une autre fonction non lineaire. par exemple une 
fonction racine, pour des calculs de distances entre des etats de neurones. Les applications peuvent aussi 
concerner d'autres dispositifs comme des g^nerateurs de fonctions, des calcuiateurs ou autres. 

Pour calculer une telle fonction. sans passer par une fonction d'approximation. on peut utiliser 
differentes manieres. 

25 On peut effectuer le calcul mathematique exact pour chaque vaieur de la grandeur independante a 
traitor, en programmant un calculateur selon les methodes connues. Une telle methode necessite d'effec- 
tuer a chaque fois les memes operations ce qui peut necesslter beaucoup de temps si le nombre de 
valeurs est eleve. 

On peut aussi prealablement stocker dans une memoire des tables precalculees. Dans ce cas, la 

30 lecture en memoire du resultat peut etre rapide. Mais pour couvrir avec un pas assez fin toutes les valeurs 
possibles de la grandeur independante, il faut alors disposer de tables de grandes capacites. Ces methodes 
de calcul presentent done des inconvenients. 

D'autre part, on peut etre conduit a identifier deux grandeurs qui sont dans la dependance I'une de 
I'autre par des couples de valeurs associant une grandeur dependante a une grandeur independante. Ainsi, 
35 dans le suivi d'un processus industriel, on peut etre conduit a mesurer par exemple un rendement R d'une 
operation en fonction de la temperature T a laquelle a ete reaiisee ladite opeation R = f(T). Un suivi du 
processus peut consister a consigner des lots de couples de mesures sur un graphique. Celui-ci peut etre 
utilise pour caracteriser le processus ou pour en deduire de nouveaux parametres de reglage de ladite 
operation. Geci fait par exemple I'objet de I'article de H. ISHIBUCHI et H. TANAKA. "Regression analysis 
40 with interval mode! by neural networks" in "IEEE International Joint Conference on Neural Networks", vol. 2.. 
18-21 Nov. 1991, SINGAPORE. II est souhaitable que ces nouveaux parametres soient representatifs du 
fondement de ladite operation et que les fluctuations de mesures inherentes a ce genre de processus 
soient ecartees. II est done souhaitable de determiner une approximation de la fonction f(.). 

Ainsi dans un cas il peut c'agir de mesures erratiques ou entachees d'erreurs que Ton desire 
45 representer par une fonction d'approximation. 

Dans un autre cas, on connait des valeurs precises mais t'utilisation a en faire ne necessite pas une 
grande precision et une fonction d'approximation suffit. 

Un des buts de I'invention est de generer une fonction d'approximation avec des moyens materiels 
reduits permettant de calculer rapidement un nombre reduit de valeurs de la grandeur dependante utiles a 
50 I'application sans avoir pour cela a determiner d'autres valeurs de la fonction d'approximation. Un but 
complementaire est de delivrer des valeurs qui peuvent etre approchees dans la limite d'une erreur 
maximale controlee. 

Ce but est atteint avec un dispositif caracterise en ce qu'il comprend ; 
- des premiers moyens : 

55 - pour determiner iterativemant au moins une fonction lineaire courante de regression en rendant 

egales, en vaieur absolue, des premieres erreurs de signes alternes mesurees entre des premieres 
valeurs de la grandeur dependante pour trois couples d'une suite desdits premiers couples, et 
respect! vement des secondes valeurs de la grandeur dependante determinees, d'apres ladite 
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fonction lineaire courante. pour les memes vaieurs de la grandeur independante. 

- pour selectionner celle des fonctions (ineaires courantes qui delivre I'approximation de tous les 
couples de ladite suite avec des erreurs minimales. 

- et pour coder, a I'aide de codes specifiques. la fonction lineaire de regression selectionnee, 

5 - et des seconds moyens pour determiner lesdits seconds couples a I'aide desdits codes specifiques. 
Ainsi avantageusement on determine une fonction lineaire de regression approchant au mreux les 
differents couples de vaieurs connus. Les resultats approches ainsi delivres forment un compromis 
satisfaisant pour de nombreuses utilisations du dispositif generateur de fonctio... 

Une fonction lineaire de regression est une fonction simplificatrice qui represente un phenomene 
10 complexe en reduisant les parametres significatifs. En representant la suite de couples de vaieurs par des 
points dans un espace a deux dimensions, la fonction lineaire de regression devient une droite de 
regression, 

Ainsi apres avoir defini la droite de regression par des codes, on peut calculer une valeur approchee de 
la grandeur dependante en tout point de la droite de regression avec des moyens reduits pour des vaieurs 
75 quelconques de la grandeur independante. 

L'invention concerne egalement une methode pour generer une fonction d'approximation, la methode 

comprenant : 

• une premiere phase : 

- pour determiner iterativement au moins une fonction lineaire courante de regression en rendant 
20 egales, en valeur absolue, des premieres erreurs. de signes alternes. mesurees entre des 

premieres vaieurs de la grandeur dependante pour trois couples d'une suite desdits premiers 
couples, et respectivement des secondes vaieurs de la grandeur dependante determinees. d'apres 
ladite fonction lineaire. pour les memes vaieurs de la grandeur Independante. 
pour selectionner celle des fonctions lineaires courantes qui delivre Tapproxlmation de tous les 
25 couples de ladite suite avec des erreurs minimales. 

- et pour coder la fonction lineaire de regression selectionnee a I'aide de codes specifiques. 

- et une seconde phase pour determiner lesdits seconds couples a I'aide desdits codes specifiques. 
Les moyens mis en oeuvre par Tinvention peuvent etre formes par un calculateur programme ou par un 
circuit dedie. Us peuvent aussi mettre en oeuvre des neurones. 

30 Un dispositif mettant en oeuvre des neurones selon l'invention peut etre utilise par un reseau de 
neurones, dont il peut notamment en constituer un sous-ensemble. En effet. pour fonctionner, le reseau de 
neurones doit disposer de moyens pour appliquer une fonction non lineaire deactivation aux potentiets de 
neurones qu'il delivre. Selon {'invention, ie dispositif muni de neurones peut calculer une approximation de 
cette fonction non-lineaire deactivation. II peut egalement calculer des distances entre des etats de neurones 

35 en calculant une approximation d'une fonction racine carree destinee a etre exploitee dans le reseau de 
neurones. 

Lorsque la taille de la suite de couples de vaieurs fournis initialement est elevee. on peut diviser la suite 
de couples en plusleurs sous-ensembles pour determiner plusieurs droltes de regression et ameliorer la 
precision de I'approximation. La fonction d'approximation de la suite de couples est alors formee par une 
40 fonction lineaire par morceaux pour laquelle une exigeance de continuite entre les morceaux peut etre ou 
non imposee. 

Certains couples de I'ensemble de couples de vaieurs peuvent avoir une influence particuliere que Ton 
peut concretiser en donnant un coefficient de ponderation specifique a chaque couple. Dans ce cas. I'erreur 
affectee a chaque couple tient compte de ce coefficient de ponderation specifique. 
45 Ces differents aspects de ('invention et d'autres encore seront apparents et elucides a partir des modes 
de realisation decrits ci-apres. 

L'invention sera mieux comprise a I'aide des figures suivantes donnees a titre d'exemples non llmitatifs 
qui representent : 

Figure 1 : un graphique montrant une representation a deux dimensions d'un ensemble de points avec 
50 une droite de regression D. 

Figure 2 : un graphique montrant un ensemble de points et des droites servant a la determination d'une 
enveloppe. 

Figure 3 : un organigramme d'une premiere varlante de mise en oeuvre de la methode a partir de 
triplets de points. 

Figure 4 : un organigramme d'une seconde variante de mise en oeuvre de la methode a partir de triplets 
de points. 

Figure 5 : une partie d'organigramme d'une troisieme variante de mise en oeuvre de la methode a partir 
de couples de points. 
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Figure 6 ; une partie d' organigram me indiquant une preselection de points appartenant a une enveloppe 
inferieure ou a une enveloppe superieure de Tensemble des points. 
Figure 7 : un organlgramme pour la determination des enveloppes inferieure et superieure. 
Figure 8 : un schema d'un dispositif selon I'lnvention. 
5 Figure 9 : un schema d*un dispositif de codage pour le calcul et le codage d'une droite de regression 
dans le cas general. 

Figure 10 : un schema d'un dispositif de codage pour le calcul et le codage d*une droite de regression 
pour la premiere variante. 

Figure 11 : un schema d'un dispositif de transcodage pour le calcul des valeurs de la grandeur 

JO dependante a partir d'un ensemble de droites de regression codees. 

Figure 12 : un schema d'une realisation neuronale des moyens de comparaison de la figure 10. 
Figure 13 : une representation montrant une fonction d'approximation formee de plusieurs droites. 
Figure 14 : deux representations montrant une determination des limites de validite de deux droites 
consecutives- 

75 Figure 15 : une representation concernant le raccordement de droites consecutives. 

L'invention concerne I'approximation d'une fonction connue uniquement a travers un certain nombre de 

points par exemple Pi... Pe (figure 1) dans une representation a -deux dimensions. Chacun de ces points 

est defini par un couple de valeurs (x. y) reliant la grandeur independante x a la grandeur dependante y. 

Par la suite, il sera question d'une paire (respectivement d'un triplet) de points, ce qui met en oeuvre deux 
20 couples (respectivement trois couples) de valeurs. On ordonne les points d'apres un ordre croissant des 

valeurs d'abscisses Xj. ce qui definit un indice i croissant avec lesdites valeurs. Une convention inverse 

peut etre faite en modifiant en consequence Texpose ci-apres. 

Selon I'invention on effectue une approximation de I'ensemble des couples (Xi, Yi). (X2. Y2)... par une 

droite de regression D ayant pour equation : 

25 

(1) D:y = p.x + Q. 

ou X et y sont des variables courantes. 

Pour cela on considere dans une premiere variante trois couples de valeurs, par exemple (X3, Y3), (X*. 

30 Y4). (Xs, Y5), et on determine une droite de regression D par equilibrage des p^reurs absolues. Une erreur 
est mesuree par la difference apparaissant. pour une abscisse x donnee. entre la valeur y du point et 
Tabscisse y mesuree sur la droite de regression. Equilibrer tes erreurs sur trois points, consiste a avoir trois 
erreurs egales en valeur absoiue avec un signe d'erreur oppose aux deux autres pour le point ayant une 
abscisse x comprise entre les abscisses x des deux autres points. Puis on examine si pour les points 

35 restants de I'ensemble. I'erreur qui les separe de la droite reste inferieure ou egale en valeur absoiue a 
Terreur prealablement determinee pour les trois points selectionnes. Ceci decoule du fait que Ton 
s'interesse a une erreur pire cas relative a I'ensemble de tous les points a prendre en consideration c'est- 
a-dire la plus grande erreur. en valeur absoiue, qui existe entre un des points et la droite de regression. Si 
toutes les erreurs sont effectivement inferieures ou egales. la droite est selectionnee pour representor les 

40 points sinon on recommence les operations avec trois nouveaux couples de valeurs pour determiner une ■ 
autre droite de reg-ession. 

II peut exister plusieurs droites de regression representant tous les points de Tensemble. Selon la 
methode on determine la droite de regression optimale qui minimise I'erreur pire cas definie prealablement. 
La figure 1 represente un exemple forme de six points Pi a Pg disposes selon une representation a 

45 deux dimensions. A des fins d'explication considerons le resultat final. On observe que la droite de 
regression D de la figure 1 est situee de sorte que les erreurs sont egales en valeur absoiue pour les points 
P3, Pa et P5. Pour les points Pi. P2. Ps les erreurs sont inferieures aux precedentes, en valeur absoiue. 
Dans un ensemble de points P^ a Ps. la methode va ainsi consister a rechercher les trois points 
particuliers. ici P3. Pa, P5, qui permettent de determiner la droite de regression optimale qui minimise 

50 i'erreur pire cas puis a coder cette droite. Dans le cas de la figure 1 representant un resultat final, si on 
trace deux droites D^ et D2. paralleles a la droite de regression D. qui passent respectivement par les 
points P3 . Ps . d'une part et Pd d'autre part, on constate que tous les points de I'ensemble sont a I'interieur 
d'un bandeau limite par les droites Di et D2 ou sur ces droites. 

La phase de determination de la droite de regression peut donner lieu a plusieurs mises en oeuvre dont 

55 seules les plus avantageuses seront decrites ci-apres. 
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Premiere vanante de la premiere phase de la methode . 

La figure 3 represente la suite des etapes a mettre en oeuvre pour determiner la droite de regression 
optimale. 

5 Parmi I'ensemble des points Pj... Pn (bloc 100). on selectionne (bloc 102) trois points quelconques Pj, 
Pj, Pk avec i < j < k. Ces trois points servent a determiner la droite de regression qui minimise I'erreur sur 
y pour ces trois points. Cette determination est effectuee de manlere analytique preferentiellement par des 
moyens programmes. On determine une droite de reg.ession D pour que trois erreurs (bloc 104) : 
EpD (Pi, D). EpD (Pj, D). Epo (Pk. D) 

10 entre la droite 0 et chacun des points verifient : 

Epo (Po D) = -Epo (Pj. D) = Epo (Pk.D) 

avec 

75 

Epo (Pi. D) = Yi - (p.Xj + q) 

et des relations analogues pour les autres erreurs. 

On determine la droite de regression D a Taide des coefficients p et q de {'equation (1 ) tels que : 

20 




25 



^ 2 {X^ - X^) 

L'erreur se rapportant a un triplet (Pi. Pj. P^) s'ecrit alors : 
35 ET{Pi,Pi,P,) = \Epo{Pi,D)\ 

Lorsque l'erreur Et a ete ainsi calculee, on examine si les autres points de I'ensemble engendrent des 
erreurs inferieures ou egales, en valeur absolue. a celles des points Pj. Pj. P^. Pour cela on selectionne un 
point (bloc 106) additionnel et on calcule la valeur absolue de l'erreur Ep„, (bloc 108) apparaissant entre 

40 la valeur de la variable y au point Pm et la droite D. 

Lorsque cette erreur Ep^ est inferieure ou egale en valeur absolue a Et (bloc 110) (repere Y). te point 
Pm additionnel est accepte et la methode se poursuit (bloc 112) avec un point additionnel suivant (bloc 
106). Si tous les points additionnels satisfont au critere | Ep^ | S Er , la droite D est accepiee et ses 
coefficients sont utilises pour coder la droite de regression optimale Dopi (bloc 114). 

45 Lorsque cette erreur Ep„, est superieure a Et (bloc 110) (repere N). le triplet de points P.. Pj, P^ 
selectionne n'est pas accepte et un autre triplet de points (bloc 116) est choisi (lien 101) dans I'ensemble 
de points (bloc 102). La methode s'acheve a I'obtention de la droite D (bloc 114) satisfaisant ce critere 
meme si tous les triplets de points n'ont pas ete examines. 

II peut apparaltre des situations ou a Tissue de I'etape 116 tous les points possibles ont ete examines 

50 et aucun des triplets n'a fourni de solution (bloc 118). Dans ce cas. il est possible de reprendre le 
deroulement de la premiere varianie en remplagant te test du bloc 1 10 par le test suivant : 

\Ep„,\S o£r 

55 ou a est un coefficient legerement superieur a i . Dans cette meme situation, il est egalement possible de 
faire appel a la deuxieme variante. 

La premiere phase peut comprendre les etapes suivantes : 
A - selection de trois couples de valeurs parmi ladite suite. 
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B - calcul de la fonction iineatre courante de regression D et determination d'une erreur de triplet 

Et=\Epd\ s'y rapportant, 
C - selection d'un couple additionnel. 

D - calcul d'une erreur additionnelle Ep„, entre le couple additionnel et ladite fonction. 
5 E - et lorsque | Ep^ | S EH110) pour le couple additionnel. la methode reprend a Tetape C avec un 
couple additionnel suivant, 
F - et lorsque | Epm \> Et pour au moins un couple additionnel, 

la methode reprend a I'etape A avec une selection d'un autre groupe de trois couples parmi ladite suite 
de couples. 

10 G ■ et lorsque \ Epn, \ S Et pour tous les couples additionnels, la fonction lineaire courante de 
regression est codee et stockee en tant que fonction lineaire d'approximation. 
On peut choisir de scruter I'ensemble des triplets en prenant un ordre croissant ou un ordre decroissant 
ou un ordre aleatoire pour effectuer cet examen. Le triplet qui sera retenu pour determiner la droite de 
regression pourra de ce fait etre detecte a un instant quelconque du deroulement de cette scrutation. II 

;5 s'ensuit que la rapidite d'obtention de la droite de regression depend de Tinstant au cours duquel le triplet 
est detecte. Sa mise en oeuvre presente un degre de complexite allant de N a ou N est ie nombre de 
points initiaux. Sa complexite est done reduite pour un petit nombre de points. Cette variante permet 
d'obtenir une realisation materielle avec une forte parallelisation. Elle est tres peu sensible a une troncature 
des valeurs et fournit un resultat exact. 
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Dans cette deuxieme variante (figure 4), on selectionne successivement des triplets de points, on 
calcule chaque fois une droite de regression, et. par recurrence, on selectionne celle qui delivre I'erreur Epo 
25 la plus grande c'est-a-dire correspondant a I'erreur pire cas pour I'ensemble de points consideres. 
Selon la seconde variante, la premiere phase comprend les etapes suivantes : 
A - selection de trois couples de valeurs parmi ladite suite, 

B - calcul d'une fonction lineaire courante de regression D et determination d'une erreur de triplet 
Et'\Epd\ s'y rapportant, 

30 C - comparaison de I'erreur Et avec une erreur optimale Eop initla'^see a une valeur striclement 
negative. 

D • et lorsque Er > Eop , mise a jour de Terreur optimale Eop en remplagant Eop par Et et mise a jour 
des codes d'une fonction lineaire optimale de regression Dop en remplagant ceux-ci par les 
codes de la fonction lineaire courante D, 
35 E - puis retour a I'etape A pour selectionner trois autres couples. 

F - et lorsque tous les triplets de couples de valeurs de la suite ont ete testes, les derniers codes de 

la fonction lineaire optimale Dop constituent les codes de la fonction lineaire d'approximation. 
A chaque examen d'un triplet, on compare I'erreur Et dudit triplet avec I'erreur optimale precedemment 
memorisee et on met a jour I'erreur optimale Eop avec la plus grande valeur de I'erreur Et determinee pour 
40 chaque triplet. On met a jour egalement les parametres de la droite optimale. Avant la mise en oeuvre de 
I'etape A. il faut initialiser la valeur Eop a une valeur faible negative, par exemple -1. 

Dans ce cas, la rapidite d'obtention de la droite de regression ne depend pas du mode de scrutation 
des triplets. Sa complexite est de degre N-. Une mise en oeuvre materielle peut beneficier de la grande 
regularite de I'algorithme mis en oeuvre. 1) est peu sensible aux troncatures des donnees et fournit une 
45 solution exacte. 

Troisieme variante de la premiere phase de la methode. 

Dans cette troisieme variante (figure 5), on selectionne d'abord une paire de points a laquelle on ajoute 
50 un point supplementaire. situe entre ces deux points, afin de former un triplet de points. Pour cela. on 
modifie les etapes A. B et C de la premiere vanante. les autres etapes restant les memes. Les etapes 
modifiees sont telles que : 

Ai • modifie I'etape A en operant une selection de deux couples de valeurs appartenant a ladite 
suite, tel qu'il existe au moins un couple additionnel intermediaire ayant une grandeur indepen- 
55 dante {X^ - X^) comprise entre les grandeurs independantes dudit couple pour constituer au 

moins un triplet de couples, 

A2 - modifie I'etape A premierement en determinant une fonction lineaire annexe qui contient les 
deux couples selectionne s et deuxiemement en determinant des secondes erreurs entre les 



6 



EP 0 624 847 Al 



grandeurs dependantes des couples intermediaires possibles et ladite fonction lineaire annexe ; 
. et lorsque ces secondes erreurs sont tootes de meme signe. selection du couple 
intermediaire fournissant la plus grande seconde erreur. en valeur absolue. pour former un 
triplet de couples de valeurs forme du couple intermediaire et des deux couples 
5 selectionnes. 

. et lorsque ces secondes erreurs sont de signes differents. reprise de la methode a Tetape 
Al. 

Bl • I'etape B est effectuee avec ledit triplet selectionne. 

Cl • modifie I'etape C en selectionnant un couple additionnel dont la grandeur independante n'est 
w pas comprise entre la grandeur independante des deux couples selectionnes. 

Lorsque Terreur | Epny | est superieure a I'erreur Et (bloc 110). la methode reprend a Tetape Al (lien 
101) avec une nouvelle selection de paire de points (bloc 102a). 

On observe que te deroulement de cette troisieme variante depend de la scrutation des valeurs et done 
des valeurs elles-memes. La complexite de la mise en oeuvre de cette variante varie entre N et N- ce qui 
15 fui confere un certain avantage par rapport aux variantes precedentes. Une forte parallelisation des moyens 
de mise en oeuvre peut etre operee mais Timplementation des moyens peut presenter un manque de 
regularite ce qui peut constituer un handicap pour une realisation integree. Cette variante est peu sensible 
aux erreurs d'arrondis des valeurs et fournit une solution exacte. 

20 Quatrieme variante de la premiere phase de ta methode. 

Elle concerne la determination de la droite de regression a partir des enveloppes. 

II est possible de reduire le nombre de triplets a examiner en determinant des enveloppes respective- 
ment ^uperieure et inferieure entourant les points extremes dans la representation bidimensionnelle de 
25 I'ensemble de points. Une enveloppe superieure ou une enveloppe inferieure est definie telle qu'en joignant 
par une droite deux points adjacents quelconques de I'enveloppe, tous les autres points soient situes d'un 
meme cote respectivement de I'enveloppe superieure ou de I'enveloppe inferieure. On determine ainsi tous 
les points appartenant a ces dites enveloppes. 

La determination de la droite de regression va consister a considerer les paires de points adjacents 
30 d'une des enveloppes auxquels on associe un point intermediaire n'appartenant pas a ladite enveloppe 
pour constituer un triplet et operer comme cela vient d'etre decrit dans le cas des paires de points de la 
troisieme variante. Si une solution optimale n'a pas ete trouvee. on considere les paires de points de Tautre 
enveloppe. 

Pour mettre en oeuvre une enveloppe. on selectionne une paire de points adjacents appartenant a 
35 I'enveloppe. On determine alors s'il existe un point dispose de telle fagon que son abscisse soit 
intermediaire entre les abscisses des points selectionnes. Lorsque ce point n'existe pas on passe a une 
autre paire de points de la meme enveloppe. Pour certaines paires. lorsqu'il apparall qu'il existe un ou 
plusieurs de ces points intermediaires. on choisit le point intermediaire le plus eloigne de la droite 
contenant les deux points de la paire pour former un triplet et pour determiner une droite de regression. 
40 Pour determiner si cette droite de regression peut etre selectionnee comme droite de regression optimale. 
la methode met en oeuvre les memes operations que celles decrites prealablement dans le cas de la 
troisieme variante. 

Pour cela on modifie la troisieme variante telle que (figure 6). prealablement a I'etape 102a (figure 5). la 
premiere phase de la methode comprend une etape (bloc lOOa) de determination d'une enveloppe 

45 inferieure et ou d'une enveloppe superieure reunissant les points les plus extremes de I'ensemble de 
points, la selection des paires de points a I'etape 102a etant faite parmi les points adjacents appartenant a 
I'une ou I'autre enveloppe. La selection de ladite paire de points est effectuee lorsqu'il existe au moins un 
point intermediaire ayant une abscisse situee entre les abscisses des points de la paire de points. S'il 
existe plusieurs points intermediaires. on forme le triplet avec le point intermediaire le plus eloigne de la 

50 droite passant par les deux points qui forment la paire de points. Si une solution n'est pas trouvee avec la 
premiere enveloppe. on poursuit le traitement avec la seconde enveloppe. 

La determination des enveloppes inferieure et superieure est effectuee selon I'organigramme de la 
figure 7. L'indice des points allant croissant avec la variable x. le premier point Pq fait ainsi partie des deux 
enveloppes. Considerons d'abord I'enveloppe inferieure. les points appartenant a I'enveloppe inferieure sont 

55 reperes par la lettre 0. Un point courant 0^ est repere par l'indice v. 

A I'etape 400 les deux premiers points de I'enveloppe sont : Go = Po et Ot = P, . Le dernier point 
courant Ov est repere : v = i pour P, . Le nombre k de points traites est comptabilise. 



EP 0 624 847 A1 



Etape 402 : un premier test est effectue pour determiner si le dernier point Pn-i a ete traite pour 
detecter la fin de la determination de I'enveloppe. 

Etape 404 : dans le cas contraire on teste si v > 1. Lorsque v < 1 on incremente v (v = v + 1) et on 
prend le point Pk comme point Qv (etape 407). On incremente Tindice k pour traiter le point P suivant 
5 (etape 409) et reprise de la methode a I'etape 402. 

Lorsque v ^ 1 . on calcule une droite passant par les points Qy - i et Qv (etape 406) et on determine le 
signe e de I'erreur sur ia variable dependante entre le point courant et cette droite (Qv -i. Qv)- Ceci a 
pour but de determiner si le point courant est au-dessus ou en dessous de la droite (Qv - i. Qv). 

Lorsque le signe « 5 o. il faut supprimer le dernier point Qv. decrementer v tel que v = v - 1 (etape 
10 410) et reprendre la methode a I'etape 404. On peut de cette maniere etre conduit a supprimer certains 
points deja acceptes lorsqu'un point suivant necessite de les supprimer. 

Lorsque le signe t est strictement positif. la methode reprend a I'etape 407 avec un point suivant. 

Get organigramme s'applique a la determination des enveloppes inferieure et superieure en inversant le 
signe de I'erreur a considerer. 
75 Pour faire comprendre les mecanismes ainsi mis en oeuvre considerons, a titre d'exemple, le cas 
simple forme par les points Pi, P2, Pa, P* de la figure 2 et determinons I'enveloppe inferieure. Le point Pi 
est te premier point de I'enveloppe d'ou Qo = Pu On prend P2 comme deuxieme point Qi . On calcule la 
droite Lo reliant Pi, P2. Le segment Pi, P2 est accepte comme premier segment de I'enveloppe basse. On 
examine si P3 est au-dessus de la droite L©, ce qui est le cas sur la figure 2. Le segment P2 P3 est accepte 
20 comme deuxieme segment de I'enveloppe basse. On calcule ensuite la droite Li reliant P2 et P3 et on 
examine si P4 est au-dessus de la droite Li . Ce n'est pas le cas sur la figure 2. Ceci entratne que le point 
P3 et la droite Li sont d'abord ecartes. De meme, P* est en dessous de la droite U, done le point P2 et la 
droite Lo sont ensuite aussi ecartes comme ne pouvant pas appartenir a I'enveloppe inferieure et les 
operations se poursuivent avec la determination de la droite L2 reliant Pi. Pa, puis I'examen des points 
25 suivants (non representes). 

La complexite de mise en oeuvre de cette variante basee sjr des enveloppes est de degre et 
depend de Tordre de scrutation des donnees. Cette complexite est moindre que celle des variantes 
precedentes et de ce fait fournit un resultat rapidement. La regularite de I'implementation est moyenne mais 
cette variante est peu sensible aux erreurs d'arrondi des valeurs et deiivre une solution exacte. 
30 Pour certaines applications. peut etre souhaitable d'accroltre la precision de ia determination de la 
fonction d'approximation dans certains domaines de la grandeur independante x et d'affecter des coeffi- 
cients de ponderation Wj aux points Pi. C'est par exemple le cas lorsque la fonction d'approximation est 
faiblement variable avec la grandeur independante x. A certains points peu vent alors etre affectes des 
coefficients de ponderation. Ceux-ci peuvent etre communs a plusieurs points ou etre individuels pour 
35 chaque point. Ces coefficients de ponderation W, sont, par la suite, consideres strictement positifs. 

Dans les cas ou des coefficients de ponderation existent, on definit une erreur Epd entre un point Pj et 
la droite de regression D tel que : 

EpoiP.D) = W,[Yi-{p.X, + q)] 

40 

ou Epo est une valeur signee. La determination de la droite de regression D pour trois points P„ Pj. P^ est 
alors modifiee en ce que les quantites p et q de I'equation (i) deviennent : 

(2) p = NUMP'DET et q = NUMQ/DET 

45 

ou les quantites NUMP, NUMQ et DET sont definies par : 

DET = Wi.Wj.(Xi-Xi) + Wi.Wk.(XK-Xj) + W..W,.(Xk-X,) 
NUMP = W,Wj.(Yj.Yi) + Wi.W,.(YK-Yj) + \N,.V\/^.{\^^r,) 
50 NUMQ = Wi.Wj.(Xi.Yi-Xi.Yj) + Wj.WK.(XK.Yj-Xj.Yk) + Wi.Wk.(Xfc.YrXi.Y,) 

Par ailleurs. I'erreur Et associee a ce triplet peut etre exprimee et calculee par ; 
avec 

55 Ej{P.,P,P,) = \Epo{P.D\ 

EpoiP.D) = -EpoiP^D) = EpD(Pk,D) 

ou 
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(3) r - 

^ DET 

5 

La premiere et la seconde vartante de ta methode de la premiere phase decrites precedemment 
peuvent etre mises en oeuvre en faisant intervenir les coefficients de ponderation ci-dessus. Cette mise en 
oeuvr*. peut etre effectuee en programmant un calculateur. 

Dans le cas ou I'on selectionne une paire de points en faisant intervenir des coefficients de ponderation 
10 W, la methode presente les amenagements suivants. 

Pour les points Pj. P|, qui forment le triplet, on definit un facteur F, qui lie les coordonnees 
respectives et les coefficients de ponderation respectifs des points Pj. P|. Pk- Au point central 1 . on associe 
le facteur F, tel que : 




20 Ce facteur F, intervlent dans ta determination du point intermediaire Pi qui sera a selectionner (1 < l< k) 
pour constituer le triplet en tenant compte des polds affectant les points. On choisit un point Pi et on 
calcule la droite de regression Di associee au triplet P,, P|. P^ et I'erreur En associee au triplet. 

Cette variante de la methode consiste a modifier uniquement I'etape A2 de la troisi^me variante' (bloc 
102b. figure 5). Cette etape determine Texistence et la valeur d*un point intermediaire servant a former un 
25 triplet. On cherche d'abord a former une droite de regression situee en dessous de P, et de Pn. Pour 
chaque point intermediaire Pi, on determine si F, = 1 et Epo (Pj. D|) < 0. Si au moins un point verifie cette 
condition, il n'existe pas de droite de regression situee en dessous de Pi et Pk. sinon on determine une 
grandeur Gmax qui est la quantite maximale parmi : 

. d'une part les quantites Epo (Pj, Dp) pour tous les points intermediaires, 
30 . d'autre part les quantites (a) suivantes. uniquement pour les points intermediaires pour lesquels F| < l 
avec 



S*il existe au moins un point intermediaire tel que F| > 1. on definit egalement une grandeur Gmin qui 
est la valeur minimale prise par (a) uniquement pour les points intermediaires tels que F| > i. On examine 
40 alors s*il existe au moins un point intermediaire tel que : 

EpD (P.. D,) S 0 

Gmax ^ cr Eti 

et Gm.n. a > Eti (a. coefficient ^ 1). SMI existe un tel point, il est choisi comme point Intermediaire pour 
former le triplet. 

45 SMI n'existe pas de point intermediaire tel que Fi > i, on teste s'il existe au moins un point 
intermediaire tel que : 
Epo (Pj. D,) ^ 0 

Gmajc S o-Eti (o : coefficient ^ 1). 

Si un tel point existe. il est choisi comme point intermediaire pour former le triplet. 
50 Si aucun triplet n*a ete constitue, on cherche a former une droite de regression situee au-dessus des 
points Pj et Pk- La meme methode est reprise en inversant le signe des erreurs Epo- 

Si aucun point P, ne peut etre selectionne, on recommence avec une autre paire Pj. Ph. 
Lorsque I'ensemble de points a traiter est trop important pour etre represente par une seule droite de 
regression, la methode code alors plusieurs droites de regression chacune etani determinee selon la 
55 methode decrite precedemment. 

La figure 13 represente un exemple dans lequel la fonction d 'approximation est constituee de plusieurs 
droites de regression. 
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Dans une premiere situation, de par la connaissance des grandeurs a traiter. on peut vouloir imposer 
des limites a chaque droite..de regression suivant la grandeur independante x. Amsi on peut vouloir 
disposer d'une droite de regression Da entre les valeurs [Xa.^J de la grandeur, borne incluse. borne 
exclue. De meme avec Db pour [X^, et Dc pour [Xc. X^f . Dans ce cas le probleme revient a determiner 
5 une droite dans un domaine linnite et a appliquer a chaque fois la naethode deja decrite. 

Mais il est possible de faire que pour chaque droite. la methode determine des limites suivant la 
grandeur independante x sans qu'elles soient fixees au depart. 

Le principe de la determination d'une limite optimale entre deux droi.es de regression adjacentes est 
represente sur la figure 14-A. Soient deux droites de regression Dl et D2 non optimales. La droite D1 est 
10 determinee a partir de Nl points et la droite 02 est determinee a partir des N2 points restants avec Nl + 
N2 = N. ou N est le nombre total de points. Les droites Dl et D2 forment une approximation avec 
respectivement des erreurs maximales El et E2. Supposons que E2 soit superieure a El et que Ton 
cherche a reduire E2. il faut alors transferer vers Dl un point appartenant a D2. Lorsque le nombre de 
points diminue, I'erreur resultante est sort inchangee soit reduite. Lorsque le nombre de points augmente. 
/5 I'erreur resultante est soit inchangee soit augmentee. En fonction du nombre n de points appartenant aux 
droites Dl et D2. les variations des erreurs El et E2 peuvent presenter les allures representees sur la figure 
14-B pour rexempie choisi. II s'ensurt que I'approximation obtenue globalement par Dl et D2 est optimale 
lorsque ces erreurs sont voisines. c'est-a-dire situees dans la zone delimitee par un trait discontinu sur la 
figure 14-B. Ces erreurs ne sont pas necessairement egales compte tenu du caractere discontinu des 
20 valeurs de la grandeur independante. 

Pour determiner ta valeur d'abscisse Xum entre les deux droites : 

- on determine la droite Dl sur un certain nombre de points et on calcule Terreur El maximale, 

- on determine la droite D2 sur les points restants et on calcule I'erreur E2 maximale, 

- on compare El et E2 et on transfere un point de la droite qui presente la plus forte erreur vers la 
25 droite qui presente la plus faible erreur, 

- on determine la valeur limite Xij^ lorsqu'il se produit une inversion dans le rapport entre les deux 
erreurs. 

Cette methode s'applique par recurrence a plusleurs droites delimitees entre elles par plusieurs limites 
d'abscisses. 

30 Les droites sont determinees a partir d'une suite discrete et limitee de grandeurs de mesure. 
Neanmoins, pour Texploitation des droites de regression, il est necessaire de definir leur domaine 
d'existence qui s'etend sur un continum de valeurs situees entre deux limites d'abscisses. Or la determina- 
tion des droites fournit une suite de droites qui ne se raccordent pas necessairement par leurs extremites. 11 
peut apparaHre utile pour certaines applications d'eviter qu'un ecart apparaisse sur la valeur de la grandeur 

35 dependante y pour des valeurs voisines (Xm^ - () et {X,im de la grandeur independante ou € est une 
tres faible valeur. II est possible de choisir que I'abscisse limite Xum appartienne exclusivement a Tune ou a 
I'autre droite. II est aussi possible de faire suivre la determination des suites de droite d'une procedure de 
raccordement de droites. 

Ceci est represente sur la figure 15. Une solution peut consister a valider la droite Dl jusqu'a I'abscisse 
40 du premier point appartenant a la droite 02 et a recalculer la droite D2 a partir de la nouvelle valeur de - 
grandeur dependante y ainsi determinee. En conservant le dernier point appartenant a D2. on determine 
ainsi une nouvelle droite D'2 representee en pointille sur la figure 15. Une procedure analogue peut etre 
appliquee pour substituer la droite D'3 a la droite 03. On obtient ainsi un ensemble forme de plusieurs 
droites de regression forme, dans Texemple decrit, par les droites 01. 0'2. 0'3. Get ensemble constitue une 
45 approximation des grandeurs de mesure en reduisant une erreur maximale entre les grandeurs de mesure 
et Tensemble de droites. Une variante plus adaptee a cette reduction consiste a faire que ce soient les 
droites correspondent aux plus fortes erreurs qui imposent leurs points limites comme points limites aux 
droites correspondant a des erreurs plus faibles. 

50 Dispositif cour une mise en oeuvre de la premiere phase 

Un schema d'un dispositif pour generer une fonction d'approximation selon {'invention est represente 
sur la figure 8. Le dispositif 5 revolt des couples (X;. Y,) associant la grandeur dependante Y, a la grandeur 
independante X,. Les couples (Xj. Y,) entrent dans les premiers moyens 10 pour determiner et coder la 
55 fonction lineaire de regression formant une approximation des couples de mesure. Les codes specifiques 
ainsi determines sont transmis (connexion 9) aux seconds moyens 17 qui determlnent des seconds couples 
(Xa. Y'a) d'apres les codes specifiques a ta fonction d'approximation. 
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On distingue d'une part la seconde variante qui determine toutes les droites de regression passant par 
chaque combinaison realisable de triplets avec la suite de premiers couples, et les autres variantes pour 
lesquelles a chaque droite de regression determinee (associee a son erreur Et). on examine si les points 
additionnels restants fournissent bien des erreurs additionnelles Inferieures ou egales a I'erreur de triplet Et. 

La figure 9 represente un dispositif adapte a la seconde variante. fl comprend : 

- une memoire 12c MEM qui stocke notamment tous les points appartenant a I'ensemble de points a 
traiter. Ces points sont representes par leurs coordonnees (x, y). et eventuellement leurs poids W ou 
leurs inverses 1/W, 

- une unite de calcul 13c COMPUT qui calcule, pour chaque triplet selectionne. la droite de regression 
adaptee a chaque triplet, c'est-a-dire les codes p. q de la droite et I'erreur Et associee au triplet. 

De plus, un controleur iic CONTR permet d'effectuer la gestion des operations et d'adresser des 
nouveaux triplets en effectuant les lectures/ecritures de la memoire 12c et le chargement de I'unite de 
calcul 13c par de nouveaux triplets. La selection de la droite de regression qui est a conserver, c*est-a-dire 
pour cette variante celle delivrant la plus grande erreur de triplet, est effectuee par I'unite de calcul 13c. 

La figure 10 correspond au cas des autres variantes pour lesquelles pour chaque droite courante de 
regression on examine si les autres points additionnels de la suite delivrent une erreur inferieure a I'erreur 
de triplet. 

Pour cela, les premiers moyens 10 comprennent une unite de calcul 13c et une unite de comparaison 
14c formant des moyens de calcul 19c. 

L'unite de calcul 13c, transmet les codes p. q de la droite de regression du triplet courant a I'unite de 
comparaison I4c COMPAR qui determine si les points additionnels de I'ensemble de points generent une 
erreur plus faible avec cette droite de regression que celle generee par les points du triplet courant.* Pour 
cela. {'unite de comparaison 1 4c effectue le test : 

I f Pm I > Er 

Si la droite de regression courante n'est pas acceptee (test positif). un autre triplet est selectionne et une 
procedure analogue est a nouveau effectuee. Si le test est negatif pour tous les points additionnels. la 
droite de regression est acceptee et ses parametres sont charges par I'unite de calcul 13c dans la memoire 
12c. 

Pour les variantes qui operent a partir d'une paire de points. I'unite de calcul 13c est programmee pour 
determiner les points intermediaires et pour former les triplets issus des paires de points, ceux-ci etant eux- 
memes eventuellement issus d'enveloppes de points. Dans ce dernier cas, I'unite de calcul I3c est 
egalement programmee pour determiner les enveloppes. L'unite de comparaison 14c effectue ensuite les 
comparaisons des erreurs se rapportant aux differents triplets. 

Pour effectuer la gestion des operations et adresser des nouveaux triplets (lorsque des triplets 
examines n'ont pas donne satisfaction), le dispositif comprend. comme precedemment. un controleur 11c 
qui gere : 

• fes lectures/ecritures de la memoire l2c. 

- le chargement de I'unite de caicul I3c par de nouveaux triplets. 

- I'appel vers I'unite de comparaison 14c de tous les points additionnels a examiner successivement. 

II est possible de mettre en oeuvre la premiere phase de la methode a I'aide de moyens de 
comparaison 14c (figure 12) ayant une structure neuronale. Considerons le cas le plus complet oij il existe 
des coefficients de ponderation associes a chaque point Pj. Lorsque ces coefficients n'existent pas il suffit 
de leur donner une valeur unite dans les explications qui vont suivre. A titre d'exemple. on peut stocker 
dans la memoire I2c (figures 9. 10) les inverses des coefficients Wj affectes individuellement aux points P,. 
L'unite de calcul 13c determine selon les equations (2) et (3), pour un triplet donne, les codes : -p, -q, Er et 
■ Et. 

La condition a tester dans l'unite I4c pour un point additionnel P„, a tester, de parametres X„. Y„, Wm, 

est : 

Y„ - pX„ - (ErWm) - q > 0 
ou 



. p.X^ ^ (ErW„) - q < 0. 
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Ces tests sont aisement realisables par une unite I4c ayant une organisation neuronale. En effet. ce 
test est constitue de fonctions lineaires. de fonctions de seuillage et de fonctions logiques aisement 
realisees dans un reseau de neurones. 

Une telle unite neuronale de comparaison est representee sur la figure 12, Elle comprend trois 

5 neurones N1, N2, N3. Les neurones Nl et N2 regoivent des donnees El, E2, E3. E4. Le neurone N3 regoit 
les sorties des neurones Nl et N2. Chacune des entrees de ces neurones est aflectee d'un coefficient 
synaptique Ci selon la technique connue mise en oeuvre dans les reseaux de neurones. Cette technique 
est par exemple decrite dans R.P. LIPr>MANN. "An introduction to computing with neural nets" IEEE ASSP 
Magazine. Apnl 1987. pp. 4 to 22. 

70 Pour effectuer les tests precedents, on programme les neurones Nl et N2 selon te tableau I : 

TABLEAU I 



15 



Neurone 


El 


Ci 


E2 


C2 


E3 


C3 


Ea 


c* 


1 




1 




- P 


I'W^ 


-Et 


1 


-q 


2 




1 




-P 




Et 


1 


-q 



20 Tous les coefficients synaptiques du neurone N3 sont egaux a 1 . 

Les neurones Nl et N2 calculent chacun un potentiel neuronal s tel que : s = E/ C,.f, 
A ce potentiel s. il faut appliquer une fonction d'activatlon A. Pour les neurones Nl et N3, la fonction 
deactivation Ai est : 

25 Ai (s) = i si s > 0 
Ai(s) = 0 si s ^ 0. 

Pour ie neurone N2, la fonction d 'activation A2 est : 

30 A2(s) = 1 si s < 0 

A2(s) = 0 si s 2: 0. 

La sortie du neurone N3 vaut 1 si la condition a tester est satisfaite et vaut 0 dans le cas inverse. 

35 On observe sur le tableau I que les donnees -p, -q, -Et. + Et sont des codes qui apparaissent comme 
coefficients synaptiques dans le schema de la figure 12. Les donnees Y„„ X„,. l/Wm sont des donnees qui 
apparaissent sur les entrees du meme schema. Lorsque tous les points additionnels opt ete testes et 
qu'une droite de regression a ete selectionnee les codes p. q. W charges dans la memoire 12c sont alors 
utilises pour mettre en oeuvre ia seconde phase de la methode (figure 1 1). 

40 L'interet presente par la realisation neuronale decrite est de pouvoir paralleliser les differentes 
operations a mettre en oeuvre selon les variantes deja decrites. Le fonctionnement d'une telle realisation 
neuronale est alors tres rapide. 

Lorsque les droites de regression ont ete determinees, la seconde phase consistant a calculer des 
seconds couples (Xa, Y'a) de valeurs des grandeurs est alors mise en oeuvre dans des moyens de 

45 decodage 17 (figure 11). Les codes des droites sont charges dans une memoire 12a qui au cours de la 
seconde phase est adressee. par un controleur 11 a, pour fournir les codes des droites de regression 
adressees. La memoire I2a, organisee en lignes par exemple, contient pour chaque droite de regression 
les parametres p, q et x,, ou xi. est la limite superieure d'abscisse pour laquelle chaque droite de regression 
est definie. 

50 La memoire 12a contient ainsi un tableau de parametres 



correspondant aux m droites de regression stockees. 



12 
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Les moyens de decodage 1 7 comprennent : 
■ le controleur 11 a. 

• la memoire l2a, 

- et une unite de decodage 13a. Celle-ci regoit una valeur de requete Xa de la grandeur independante 
5 pour laquelle on desire obtenir un resultat Y'a issu de la fonction d'approximation. Pour cela le 

controleur ila presents successivement certaines llgnes de codes de la memoire 12a et, pour chaque 
ligne de codes, i'untte de decodage 13a teste le code pour determiner si la valeur Xa de la 
grandeur d'entree peut etre traitee par la droite de regression codee par ladite ligne. Si cela n'esi pas 
le cas, I'unite l3a requiert du controleur iia une autre ligne de codes. Lorsque cette condition est 
10 verifiee, I'unite I3a effectue le calcul : 

Y'a = Pn-XA + qn 

si la ligne d'ordre n est la ligne adequate. 
15 Comme cela vient d'etre decrit. les premiers moyens 10 peuvent constituer un dispositif de codage et 
les seconds moyens 17 peuvent constituer un dispositif de decodage, les deux dispositi^s etant separes. lis 
peuvent en effet etre distants Tun de I'autre, les operations de codage etant separees des operations de 
decodage. Mais les deux moyens 10 et 17 peuvent aussi former un dispositif de codage/decodage unique. 
Dans ce cas certains organes peuvent etre affectes successivement a plusieurs operations. Notamment : 
20 • il peut s'agir d'une memoire 12 unique qui peut stocker les parametres des couples de vaieurs et les 
codes des droites de regression determinees. Les unites I2a et 12c forment alors la memoire 12, 

- il peut s'agir d*un controleur 1 1 unique qui regroupe les fonctions des controleurs 11c et 1 la, 

- I'unite de calcul I3c et I'unite de decodage I3a peuvent egalement former une seule unite 13. 

II est egalement possible de realiser les moyens de decodage 17 (figure 11) selon une organisation a 

25 structure parallele. Une solution consiste alors a comparer, en parall^le, la grandeur de requete Xa h tous 
les codes xl- On obtient un ensemble de signaux Tl representant les resultats de ces tests de comparaison. 
On realise ensuite. en parallele, pour toutes les vaieurs de L. un test qui combine deux a deux les signaux 
Tl et Tl..i et qui indique si la grandeur de requete Xa est dans Tintervalle d'indice L Ce test est verifle pour 
une unique valeur de L. notee n. On accede alors aux codes Pn et q„ correspondants, et on en deduit la 

30 grandeur Y'a de la meme maniere que precedemment. 

Les premiers moyens 10 (de codage) associes aux seconds moyens 17 (de decodage) peuvent etre 
utilises pour determiner la valeur d'une fonction d'approximation par au moins une droite de regression. 
Ces determinations peuvent etre effectuees pour des vaieurs queiconques de la grandeur independante (a 
rinterieur de bornes predefinies definissant I'etendue de Taction de chaque droite de regression). Cette 

35 methode evite de stocker inutilement des tables de vaieurs. pour lesquelles toutes les vaieurs ne seront pas 
exploitees. Selon la methode on ne determine que les seules vaieurs necessaires a I'application. La 
methode selon ['invention presente I'interet de ne calculer que les vaieurs necessaires. Le dispositif de 
I'invention peut etre utilise en combinaison avec un processeur neuronal pour calculer une approximation 
d'une fonction non lineaire par exemple une fonction sigmoVde. Un processeur neuronal calcule des 

40 potentiels neuronaux Issus de sommes de produits de coefficients synaptiques par des etats de neurones. 
La technique est developpee dans le document de R.P. LIPPMANN deja cite. Chaque potentiel neuronal 
dolt etre soumis a Taction d'une fonction non lineaire. Pour cela. chaque potentiel neuronal constitue la 
grandeur independante decrite precedemment et le calcul de la fonction th(x) par exemple est effectue pour 
les seules vaieurs de requetes utiles. La methode fournit des vaieurs approximees, mals cette approxima- 

45 tion peut etre obtenue avec une precision qui peut etre accrue selon le nombre de droites de regression 
calcuiees et eventuellement les coefficients de ponderation W. La precision reste neanmolns tributaire de la 
precision des coordonnees des couples de depart. 

Une telle methode est particulierement interessante pour les calculs de fonctions connues (telles que 
des fonctions mathematiques) ou de fonctions inconnues expllcltement. par exemple une fonction represen- 

50 tee par des points de mesures. que Ton desire slmplifier par une fonction lineaire de regression. 

Dans sa version neuronale, Tinventlon est interessante dans des applications neuronales car elle 
apporte non seulement un traitement homogene mals aussi une grande compacite de Tarchitecture 
materielle necessalre. 

II est possible de combiner les variantes de la methode. dans leurs mises en oeuvre materielles pour 
55 combiner leurs avantages. Ainsi par exemple, les combinaisons sulvantes sont avantageuses : 
Dans le cas ou les poids existent, deux combinaisons apparalssent Interessantes : 
Dans la premiere combinaison, on utilise tout d'abord la premiere variante basee sur des triplets de 
points. Si une solution n'est pas trouvee. on enchaine avec la deuxieme variante basee aussi sur des 
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triplets de points. Ceci permet d'obtenir les avantages suivants : 
Queues que soient les donnees, on obtient une solution. 

Dans le cas ou on exige qu'on obtienne toujours une solution exacte, cette combinaison permet 
d'obtenir ia rapidite nnaximale. En effet. dans le cas general la solution est fournie directennent par 
5 la premiere variante appelee, qui est celle fournissant le plus rapidement une solution exacte dans 

le cas oij les poids existent. La deuxieme variante est beaucoup plus lente mais elle permet de 
garantir que la combinaison fournit toujours une solution. 
Dans la deuxieme combinaison, on utilise tout d'abord la variante basee sur des paires de points. Si 
elle ne trouve pas une solution (qui peut etre approchee). on utilise la deuxieme variante basee sur 
JO des triplets de points. Ceci permet d'obtenir les avantages suivants : 

Quelles que soient les donnees. on obtient une solution (eventuellement approchee). 
. Cette combinaison est en moyenne plus rapide que la precedente. 
. En contrepartie. elle ne fournit generalement qu'une solution approchee. 
Dans le cas ou les poids n'existent pas. une combinaison interessante consiste a utiliser tout d'abord la 
?5 variante basee sur I'enveloppe des points. Si elle ne trouve pas de solution, on utilise la deuxieme variante 
basee sur des triplets de points. Ceci permet d'obtenir les avantages suivants : 
Quelles que soient les donnees, on obtient une solution, 
Cette combinaison permet d'obtenir la meilleure rapidite moyenne. 
Lorsque la fonction lineaire de regression a ete determinee selon la methode decrite. il est possible de 
20 determiner non seulement des grandeurs Y'a ressortissant a ladite fonction mais egalement des grandeurs 
ressortissant a la derivee de ladite fonction. En effet. pour une valeur de requete Xa de la grandeur 
ind^pendante, il suffit d'extraire des moyens de memorisation ta valeur p ressortissant a la valeur Xa pour 
connattre la valeur de cette derivee correspondante. 

25 Revendications 

1. Dispositif (5) pour generer une fonction d'approximation fondee sur des premiers couples ((Xi. Yi) a 
(Xs. Yg)) de valours associant une grandeur dependante (Yi a Ye) a une grandeur independante (Xi a 
Xe), et pour determiner des seconds couples (Xa, Y'a) de valeurs desdites grandeurs d'apres ladite 
30 fonction d'approximation caracterise en ce que le dispositif comprend : 

- des premiers moyens (10) : 

- pour determiner iterativemant au moins une fonction lineaire courante de regression en 
rendant egales. en valeur absolue. des premieres erreurs. de signes alternes, mesurees entre 
des premieres valeurs (Y3. Ya. Y5) de la grandeur dependante pour trois couples (X3, Y3) (X4. 

35 " Yft) (X5, Y5) d'une suite desdits premiers couples, et respectivement des secondes valeurs 

(Y'3, Y'a, Y*=) de la grandeur dependante determinees, d'apres ladite fonction lineaire courante. 
pour les memos valeurs (X3, X*. X5) de la grandeur independante. 

- pour selecttonner ceile des fonctions lineaires courantes qui delivre I'approximation de tous les 
couples de ladite suite avec des erreurs minimales. 

40 ' et pour coder, a I'aide de codes specifiques (p, q), la fonction lineaire de regression • 

selectionnee, 

- et des seconds moyens (17) pour determiner lesdits seconds couples (Xa. Y'a) a I'aide desdits 
codes specifiques. 

Dispositif selon la revendication 1 caracterise en ce que les premiers moyens (10) comprennent ; 

- des moyens (12c) pour memoriser les grandeurs des premiers couples et lesdits codes specifi- 
ques. 

- des moyens (13c) (19c) pour calculer successivement les codes specifiques (p, q) de fonctions 
lineaires courantes de regression et leurs premieres erreurs. egalisees en valeur absolue. 
respectives (Ey) et pour selectionner au moins une fonction lineaire courante de regression. 

- et des moyens (lie) pour controler les moyens (12c) pour memoriser et les moyens (13c) pour 
calculer. 

3. Dispositif (5) selon les revendications 1 ou 2 caracterise en ce que pour selectionner ladite fonction 
55 lineaire, les premiers moyens comprennent des moyens (I4c) : 

- pour determiner pareillement des erreurs additionnelles pour les autres couples (Xi. Yi) (X2. Y2) 
(Xg, Ye) de valeurs de ladite suite. 

- pour comparer les erreurs additionnelles avec lesdites premieres erreurs egalisees (Et). 



45 2. 
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- et pour selectionner la fonction lineaire courante qui fournit des erreurs additionnelles inferieures 
ou egales en vateur absolue auxdites premieres erreurs egalisees (Et). 

Dispositif selon la revendication 2 caracterise en ce que, pour selectionner la fonction lineaire de 
regression, les moyens pour calcuier (I3c) selectionnent celle des fonctions lineaires courantes de 
regression qui delivre des premieres erreurs egalisees (Et) maximales. 

Dispositif selon une des revendicatrons 1 a 4 caracterise en ce que les seconds moyens (17) 
comprennent : 

- des moyens (I2a) pour memoriser lesdits codes specifiques d'au moins une fonction lineaire de 
regression. 

• des moyens (13a) pour selectionner, sur requete d'une valeur de grandeur independante (Xa), 
une fonction lineaire de regression ressortissant a ladite valeur. et pour calcuter, a partir des 
codes specifiques reQus des moyens pour memoriser. une valeur de grandeur dependante (Y'a) 
associee a ta valeur de grandeur independante (Xa). 

• et des moyens (11a) pour adresser aux moyens pour calcuier, et tjr leur demande. les codes 
specifiques des fonctions lineaires de regression. 

Dans un dispositif pour generer une fonction d'approximation fondee sur des premiers couples ((Xi . Yi ) 
- (Xs. Ys)) de valours associant une grandeur dependante (Yi - Ye) a une grandeur independante (Xi - 
Xe). methode pour generer ladite fonction d'approximation a partir des premiers couples (X:. Yt) (Xg. 
Ye) de valours, et pour determiner des seconds couples (Xa, Y'a) de valours desdites grandeurs 
d'apres ladite fonction d'approximation, caracterisee en ce que, la methode comprend : 
une premiere phase : 

- pour determiner iterativement au moins une fonction lineaire courante de regression en 
rendant egales, en valeur absolue. des premieres erreurs Epo, de signes atternes. mesurees 
entre des premieres valours (Y3, Y*. Y5) de la grandeur dependante pour trois couples (X3, 
Y3). (X*, Y4) (Xs, Ys) d'une suite desdits premiers couples, et respectivement des secondes 
valeurs (Y'3, V'l, Y's) de la grandeur dependante determinee, d'apres ladite fonction lineaire, 
pour les memos valeurs (X3, X*. Xs) de la grandeur independante, 

- pour selectionner celle des fonctions lineaires courantes qui delivre I'approximation de tous les 
couples de ladite suite avec des erreurs minimales. 

- et pour coder la fonction lineaire de regression selectionnee a I'aide de codes specifiques, 

• et une seconde phase pour determiner lesdits seconds couples (Xa. Y'a) a I'aide desdits codes 
specifiques. 

Methode selon la revendication 6 caracterisee en ce que la premiere phase comprend les etapes 
suivanles : 

A - selection (102) de trois couples (P3. P4. P5) de valeurs parmi ladite suite. 

B - calcul (104) de la fonction lineaire courante de regression D et deternnination d'une erreur de 

triplet Er = |Epd| s'y rapportant. 
C - selection (106) d'un couple additionnel (Pi), (P2). (P6). 

D - calcul (108) d'une erreur additionnelle Ep„, entre le couple additionnel et ladite fonction. 

E - et lorsque | Ep„ | < £7(1 10) pour le couple additionnel. la methode reprend a I'etape C avec 

un couple additionnel suivant. 
F - et lorsque | Epm \ > Et pour au moins un couple additionnel (110). la methode reprend a 

I'etape A avec une selection d'un autre groupe de trois couples parmi ladite suite de 

couples, 

G - et lorsque | fp^ | < Et pour tous les couples additionnels (112), la fonction lineaire courante 
de regression est codee et stockee (114) en tant que fonction lineaire d'approximation. 

Methode selon la revendication 6 caracterisee en ce que la premiere phase comprend les etapes 
suivantes : 

A - selection (102) de trois couples (P3. P4. P5) de valeurs parmi ladite suite. 

B - calcul (104) d'une fonction lineaire courante de regression D et determination d'une erreur de 

triplet Er=|Epo| s'y rapportant. 
C- comparaison (210) de I'erreur Et avec une erreur optimale Eop initialisee a une valeur 

strictemeni negative. 
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D - et lorsque Et > £op - mise a jour (1 14) de I'erreur oplimale Eop en remplagant Eop par Et et 

mise a jour des codes d'une fonction lineaire optimale de regression Dop en remplagant 

ceux-ci par les codes de la fonction lineaire courante (D). 
E - puis retour (1 16) a I'etape A pour selectionner trois autres couples (102). 
F - et lorsque tous les triplets de couples de valeurs de la suite ont ete testes (116). les derniers 

codes de la fonction lineaire optimale Dop constituent les codes de la fonction lineaire 

d'approximation (114). 

9. rs/Iethode selon la revendication 8 caracterisee en ce que la premiere phase comprend les etapes 
modifiees suivantes : 

Al - modifie I'etape A en operant une selection (102a) de deux couples de valeurs appartenant a 
ladite suite, tel quNI existe au moins un couple additionnet intermediaire ayant une grandeur 
independante (Xi - Xg) comprise entre les grandeurs 
independantes dudit couple pour constituer au moins un triplet de couples de valeurs. 

A2- modifie I'etape A premierement en determinant (102b) une fonction lineaire annexe qui 
contient les deux couples selectionnes et deuxiemement en determinant des secondes 
erreurs entre les grandeurs dependantes des couples intermediaires possibles et ladite 
fonction lineaire annexe : 

et lorsque ces secondes erreurs sont toutes de meme signe, selection du couple 
intermediaire fournissant la plus grande seconde erreur. en valeur absolue, pour 
former un triplet de couples de valeurs forme du couple intermediaire et des deux 
couples selectionnes, 

et lorsque ces secondes erreurs sont de signes differents. reprise de la methode a 
I'etape Al , 

Bl '- I'etape B (104) est effectuee avec ledit triplet selectionne, 

CI - modifie I'etape C (106) en selectionnant (106a) un couple additionnel dont la grandeur 
independante n'est pas comprise entre la grandeur independante des deux couples selec- 
tionnes. 

10. Methode selon la revendication 9 modifiee en ce que. prealablement a I'etape Al, lesdits couples de 
valeurs etant representes par des points (Pi - P6) dans un espace a deux dimensions, la premiere 
phase de la methode comprend une etape (100a) de determination d'une enveloppe inferieure eVou 
d'une enveloppe superieure reunissant les points les plus extremes de Tensemble de points, la 
selection des deux couples de valeurs a I'etape Ai etant faite parmi les points adjacents appartenant a 
Tune ou {'autre enveloppe. 

11. Methode selon une des revendications 6 a 10 caracterisee en ce qu'a chaque valeur. prise dans la 
suite de valeurs de la grandeur independante. est associee un coefficient specifique de ponderation qui 
pondere lesdites premieres erreurs. 

12. Methode selon une des revendications 6 a 11 caracterisee en ce qu'elle genere une fonction 
d'approximation comprenant une suite de fonctions lineaires de regression raccordees entre elles. 
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